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Although it has been believed that the models with generalised Chaplygin gas do not 
JH , contain singularities, in a previous work we have studied how a big freeze could take 

place in some kinds of phantom generalised Chaplygin gas. In the present work, we 
study some types of generalised Chaplygin gas in order to show how different sorts of 
singularities could appears in such models, in the future or in the past. We point out 
that: (i) singularities may not be originated from the phantom nature of the fluid, and 
(ii) if initially the tension of the brane in a brane-world Chaplygin model is large enough 
then an infrared cut off appears in the past. 

Keywords: Dark energy, Future singularities 

1. Introduction 

In the past few years the observational data have forcefully shown that the current 
evolution of the universe is accelerating, so prompting a large number of the- 
oretical models aiming at explaining these data, including quintessence scenarios, 

[9] 

modifications of gravity, models which contain extra dimensions, etc ^ One of the 
most popular among these models is the generalised Chaplygin gas (GCG) which 
was initially introduced motivated by its connection with string theory as, fur- 
thermore it can admit a supersymmetric generalisation - . In addition, GCG could 
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explain a smooth transition from a dark-matter-dominated phase to the current 
epoch, as it interpolate in a natural way between dust at small scale factors and 
a cosmological constant on large scale factors. The possibility that the stuff which 
causes the accelerating expansion of the universe is a phantom fluid, i. e., a fluid 
with the equation of state parameter w < — 1 (where p = wp with p > 0) , has led 
to the consideration of the phantom GCG (PGCG) models. Although by itself a 
phantom fluid could produce the known big rip singularity i n th e future SI it has 
been shown that PGCG prevents that doomsday take place ^El. Even so, PGCG 
models are not free from other future singularities; in fact, some of such models have 
been recently seen to show the so-called big freeze singularity which is characterised 
by a finite value of the scale factor at a given finite time, On the other hand, 
a dual phantom fluid can also be defined on a Randall-Sundrum type 1 scenario 
^. Such a fluid is characterised by an equation of state parameter w < — 1 when 
the null energy condition is fulfilled. It appears therefore of interest to study the 
dual PGCG (DPGCG) ^ as well. Also considered in the present paper will be the 
different kinds of GCG and their characteristics, emphasising those models which 
show past or future singularities. It is obtained that most of the studied scenarios 
appear plagued by distinct kinds of singularities both in the future and in the past. 

This paper can be outlined as follows. In section II we analyse the different 
PGCG models and show their different behaviours in the future and the past, clas- 
sifying the involved singularities. DPGCG models are considered in section III, 
checking that whereas the brane evolution describes a universe whose future evolu- 
tion is similar to those considered in section II; in the past such a evolution could 
drastically differ, depending on the value of the brane tension. We see in section 
IV that even in the absence of a phantom fluid, past and future singularities could 
appear in the case of a GCG satisfying null, weak and strong energy conditions. In 
section V we summarise and conclude, adding some further comments. 



2. PGCG 

The PGCG model was originally introduced in Refs. 1^121. ft satisfies the same equa- 
tion of state as GCG^l, i. e. 

where A is a positive constant and a is a parameter. If a = 1, Eq. ([I]) corresponds to 
the equation of state of a Chaplygin gas. The conservation of the energy momentum 
tensor then implies 

p= { A+ ^y +a > (2) 

with B a constant parameter. 

In order to consider a phantom GCG (PGCG) two conditions have to be satisfied: 
(i) p > and (ii) p + p < 0. One can notice that there are four possibilities: 
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(i) A > 0, B < and 1 + a > 0. 

(ii) A > 0, B < and 1 + a < 0. 

(iii) A < 0, B > 0, 1 + a > and (1 + a) -1 = 2n, with n some positive integer 
number. 

(iv) A < 0, B > 0, I + a < and (1 + a) -1 = 2n, with n some negative integer 
number. 

carrying out these conditions for a fluid described by Eqs. |T]) and @. Each of 
these possibilities corresponds to a different type of PGCG and can be considered 
as a different cosmological model. Condition (ii) for a universe whose energy density 
budget is dominated by phantom energy results in a super-inflationary expansion, in 
fact a super-accelerating phase. The different sorts of PGCG fulfil the conservation 
of the energy momentum tensor. 



2.1. PGCG I 

For this model A > 0, B < and 1 + a > 0. This fluid has been studied in great 
detail in Ref. El and we review here the characteristics of a universe filled with it, 
both in the past and in the future. 

From Eq. ([2]), we can realise that the cosmic scale factor takes values in the 
interval a m i n < a < oo which corresponds to < p < A x ^ 1+a \ where 



3(l + <») 

(3) 



Then, in this model the scale factor grows up from a non- vanishing value to infinity, 
whereas the energy density increases as the universe expands until it reaches its 

2 

asymptotic value A 1 ^ 1+a K Throughout this paper, we use units such that -jf_= 1 
and H = c = 1, where k\ is the gravitational constant. As it is shown in Ref. El we 
can integrate analytically the Friedmann equation (where here and in what follows, 
we choose the fiat chart to describe the geometry of a homogeneous and isotropic 
universe) to note that the cosmic time required by the scale factor to be infinitely 
large is infinite. Therefore, the u niver se would avoid a big rip singularity SI, even 
if it is filled with phantom energy I^El Hence, the universe would expand infinitely 
in the future cosmic time, reaching an infinite size with a finite asymptotic energy 
density; i.e the universe would be asymptotically de Sitter El. 

The time elapsed since the origin of the universe to a given scale factor close to 
its initial value can be easily estimated by analysing the behaviour of the model in 
the neighbourhood of a m i n . For ao (at a given cosmic time to) close to but still larger 
than a m j n , the scale factor dependence on the cosmic time, t, can be approximated 
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by 



«o 



- 1 



1 + 2a 
2(1 + a) 



(3A(l+a))TOT (t- to ) 



2(1 + °) 
l + 2a 



for a 7^ —1/2 and by 



1 



a 



1 exp 



;A(t-*o) 



(4) 



(5) 



for a = —1/2. Then, if — 1 < a < —1/2, the universe ought to have been born 
an infinite time ago, avoiding the standard big bang singularitjfl In fact, both the 
energy density and the pressure vanish at the infinite past. On the other hand, if 
— 1/2 < a we can also estimate the time interval since the universe was born until 
it reaches a size with scale factor ao, not much larger than a m i n , to be 

_J 1 + 2q 

2(1 + a) ' 



to — t n 



1 



«0 



- 1 



(6) 



3A(l + a)J l + 2a 

The energy density, but not necessarily the pressure, is again well defined at t min 
(where it vanishes). The question becomes now appropriate, does the universe face a 
singularity at this minimum scale factor? This is the question we tackle next. Using 
the Friedmann equation, the conservation law of a perfect fluid and Eq. ([2]), the 
time derivative of the Hubble parameter, H, and the scalar curvature for a GCG 
can be written as 



p), 



(J) 



D 



B 



Consequently, if < a the pressure, and both these two quantities blow up, i.e. there 
is a past sudden singularity a big brake singularity ^ or a Type II singularity 
in the notation of Ref. with the peculiarity that the energy density vanishes 
rather than acquiring a positive non- vanishing value. But if —1/2 < a < 0, even 
though R, H, H and the pressure are all well defined at a m i n , higher derivatives of 
the Hubble parameter could still diverge. In fact, the n-th derivative of the Hubble 
parameter respect to the cosmic time has the following structure 

— H OC + Clf} -(.n-l)(±+ a ) + l + higher pQwerg in ft ( 9 ) 

where C\ is a constant. From this expression, we conclude that the n-th derivative 
diverges at a m i n if 1/(1 + 2a) < n, unless a can be written as a — l/(2p) — 1/2, 



a Although there is no standard big bang singularity as the scalar curvature is well denned, the 
model is singular in curvature for a non cosmological observer Li2J We thank L. Fernandez-Jambrina 
for pointing out this to us. 
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where p is a positive integer. For a given l/(2p) — 1/2), the n-th derivative of 
H where n = 1 + E(l/(2a + 1)) reaches very large values at a m i n . In the previous 
expression E denotes the integer valu e fu nction. This corresponds to a type IV 
singularity in the terminology of Ref. ^21. In addition, this singularity is quieter 
than a sudden singularity as the energy density and pressure are well defined at the 
singularity, as H and its first time derivative also are, but higher derivatives of the 
Hubble parameter can blow up. In the particular case a = l/(2p) — 1/2, where p is 
an integer, all the derivatives of the Hubble parameter respect to the cosmic time 
are always well defined, however. The reason is the following: the 1, .., (p — l)-th 
derivatives are well behaved at i m i n (see Eq. (O), the next derivative, which reads 

d p p±i 

-j-^H oc C 2 + P 2p + higher powers in p, (10) 

with C*2 a constant, which is also well defined at the minimum size of the universe. 
Then 

d p+1 

+1 H cx p ' + higher powers in p. (11) 

Therefore, the lower order of p in the (p+l)-th derivative is similar to the order of p 
in H and so well defined, too. Thus, we are back to our starting point. We can then 
conclude that all the derivatives of H are well defined at i nnn if a = l/(2p) — 1/2 
being p a positive integer. 

In summary, if— 1 < a < —1/2, then there is no singularity at a finite time in the 
past. If —1/2 < a and a l/(2p) — 1/2, where p is an integer, there is a singularity at 
amin which can be of two different types depending on (i) < a or (ii) — l/2<a<0. 
In the former case, i.e. < a, it corresponds to a sudden past singularity with 
a vanishing energy density, while in the latter case, —1/2 < a < 0, it corresponds 
to a divergence in higher derivativ es o f the Hubble parameter (in fact, a type IV 
singularity^ in the notation of Ref. 021). 



2.2. PGCG II 

Now, we consider a PG CG like in the previous case; i.e. A > and B < 0, but 
where 1 + a < l ^^ l Here the scale factor is bounded < a < a max and the 
energy density is larger than ^4 1 /( 1 + Q ), where 

B ^ofer 



It is a singularity as far as the derivatives of the Hubble parameter diverge, even though the 
scalar curvature and Riemann tensor do not. 
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The cosmic time scales with the scale factor as 

2 



tr 



t : 



.^4-2(1 + 0) 



3(1 + 2a)' 

1 + 2a 



1 



-3(l+a) 



l + 2a 
2{l + c) 



x F 1, 



3 + 4a 



2(1 + a)' 2(1 + a) 



1 



-3(l+a)^ 



(13) 



where F(6, c; d; e) is a hypergeometric function El i n the previous expression t max 
corresponds to the cosmic time when a universe filled by this fluid would hit a 
big freeze singularity^ i.e. a singularity at a finite scale factor, a max , and cosmic 
time, t max , and where both the energy density and pressure blow up, as the Hubble 
parameter and its cosmic derivative^ (see Fig.Q]) do as well. It can be easily checked 
that the cosmic time since the universe has a given size (at a given t) till it hits the 
big freeze is finitcQ . This result confirms our previous approximation close to cimax 

m 



a ~ a n 



1 + 2a 



2(1 + a-) 



2(1 + °) 
l + 2a 



A— 13(1 + a)\^(t r , 



2(1 + Q ) 
t) 1 + 2 = 



(14) 



It can also be noticed that when the scale factor reaches its maximum finite size, the 
scalar curvature blows up (see Eq. ©). This sin gularity can be classified as Type III 
in the notation of Ref. H^l (see also Refs. I ^^ l ^^ l). In this model the universe expands 
until it reaches a maximum size, at which moment it freezes with an infinite energy 
density; for this reason this doomsday has been named big freeze singularity 1^1. 

Finally, let us briefly point out that a universe filled with this kind of PGCG is 
asymptotically de Sitter in the past. The scale factor in the regime of small scale 
factor can be approximated by 



ao exp 



(Aw+")(t-t )) 



(15) 



where ao is a non-vanishing small scale factor. Therefore, the universe starts its 
evolution at an infinite past (a — * as t — t Q —> 00) where there is no big bang 
singularity. 



2.3. PGCG III and IV 

These PGCG correspond to having A < and B > in Eqs. ([I])-©. Hence, in order 
to satisfy the phantom energy conditions; i.e. p > and p + p < 0, the parameter a 
must be quantised so that (1 + a)~ x — 2n, with n an integer. There are two kinds 
of phantom fluids depending on n; i.e. (1 + a)^ 1 , positive or negative. We will refer 



c A hypergeometric series F(fe, c; d; e), also called a hypergeometric function, converges at any value 
e such that |e| < 1, whenever b + c — d < 0. However, if0<6 + c — d < I the series does not 
con verge at e = 1. In addition, if 1 < b + c — d, the hypergeometric function blows up at |e| = 1 
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Fig. 1. A schematic plot of a future big freeze singularity. This singularity takes places at a finite 
future cosmic time where not only the energy density but also the pressure blows up. Therefore, 
both the Hubble rate and its cosmic time diverge as well at the big freeze event. 



to the former fluid (n > 0) as the phantom GCG III and to the latter (n < 0) as 
the phantom GCG IV. 

A universe filled with a type III phantom GCG would have a minimum scale 
factor (a m i n < a, where a m ; n is given by Eq. ((3])) at an infinite past which can be 
noticed by using the conveniently modified Eqs. @i and when A is replaced by 
\A\. As it can be realised, the advantage of having a quantised a is to remove the 
possible different past singularities (at a finite cosmic time) present in a phantom 
GCG I model. After all, this advantage is not that surprising as having n > 
implies —1 < a < —1/2, precisely a situation without past singularities (at a finite 
cosmic time) for the phantom GCG I . Finally, we note that such a universe would 
be asymptotically de Sitter in the future. 

On the other hand, a universe filled with a type IV phantom GCG would have 
a maximum scale factor (a < a max , where a max is given by Eq. (jT3J)) . Then, the 
universe would hit a big freeze singularity at a finite future cosmic time (Eq. (|14[) 
after replacing A for \A\ applies in this case), where the energy density and pressure 
would diverge (see Fig. QJ. Therefore, unlike the previous case, the quantisation of 
the parameter a is useless to remove the big freeze singularity. Let us point out that 
such a universe would be asymptotically de Sitter in the past and born an infinite 
time ago. 

A summary of the results found on this section are shown in table I. 



3. DPGCG 

In this section, we thoroughly analyse the cosmology of a dual phantom GCG 
(DPGCG) Such a dual fluid is defined as a perfect fluid which satisfies the 
equation of state JT]) with w < — 1 and fulfils the null energy condition' ' ' ' 10 '. This 
definition of dual has not to be confused with the one used in Ref. 

The key point about the dual of a phantom energy setup is that it may supply 
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Table 1. This tabic summarises the asymptotic behaviour of a universe filled with each type of a phantom GCG. The past asymptotic 
behaviours labelled by (1) and (3) correspond to —1 < a < — 1/2 and < a, respectively. The past asymptotic behaviour labelled by 
(2) corresponds to— l/2<a<0, where a cannot be expressed as a = l/(2p) — 1/2, with p a positive integer. If— 1/2<q<0 and a 
can be expressed as a = l/(2p) — 1/2, with p a positive integer, there is no past singularity and the universe is born at a finite past. 



A,B 


1 + a 


a 


P 


Past 


Future 


A > 
B < 


positive 


Gmin < a < OO 


< P < A 1 /(i+«) 


(1) infinite past 


asymptotically de Sitter 


(2) type IV singularity 


(3) sudden singularity 


negative 


< a < Clmax 


AVU+«J < p < oo 


asymptotically de Sitter/infinite past 


big freeze singularity 


A < 
B > 


(2n)- L > 


Gmin < a < CO 


< P < \A\ L ^ 1+a > 


infinite past 


asymptotically de Sitter 


(2n)- 1 < 


< a < a max 


| A |i/ti+a) < p < oo 


asymptotically dc Sitter/infinite past 


big freeze singularity 



an alternative for dark energy mo dels in theories with modified Friedmann equation 
like the brane-world scenario des pite its energy density is negativ^. We next 
review how this can be the case in a Randall-Sundrum model with a unique 
brane (RSI) ^1 The modified Friedmann equation on the brane reads 



(16) 



where A is the positive brane tension. Therefore, the square of the Hubble parameter 
in RSI scenario is well defined even for a negative energy density p as long as 
p < — 2A; i.e. as long as the effective energy density of the brane is positive. In this 
setup, the brane is super-inflating; i.e. H > where 



* = -§(p + ,)(i + £)>o. 



(17) 



Consequently, a brane filled with the dual of a phantom energy is always acceler- 
ating; i.e. a > 0, in particular if it is filled with DPGCG. Moreover, the brane is 
super-inflating although the dual of the phantom energy satisfies the null energy 
condition; i. e. p + p > 0, because p is negative and smaller than — 2A. We remind 
that here the effective energy density of the brane is positive and is well behaved. 
We will analyse three different types of DPGCG characterised bjEI 



(i) A < 0, B > 0, 1 + a > and (1 



2n + 1, with n a positive integer, 



d A negative dark energy component has been also used in Ref. but in different set up. 
c The families of duals of the phantom GCG is much larger and can be characterised by: 

1 + 2(m + n) 



1 + a 



1 + 2n 



-, m,n e Z, A < 0,B > 0, 

l + 2n 



B 



3[l + 2(n 



a 



l + 2(m + n) 



< 0, with - \A\ + 



3[l + 2(m + n)] 

a !+ 2 ™ 



< 0, 



P=\A\[ \A\ 



3[l + 2(m + n)] 

a !+ 2 " 



> 0. 



(18) 

(19) 
(20) 



Here we restrict to the case m = —n. These set of solutions fulfil the conservation of the energy 
momentum tensor. 
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(ii) A < 0, B > 0, 1 + a < 0, (1 + a) -1 = 2n+ 1, with n a negative integer, and 
|A|Vd+o) < 2A, 

(iii) A < 0, B > 0, 1 + a < 0, (1 + a) -1 = 2n + 1, with n a negative integer, and 
2A < |A|V(i+«), 

which we will refer to as DPGCG I, II and III, respectively. As it can be noticed 
a commonness of DPGCGs is that the parameter a is quantised; more precisely, 
1 + a = 1/(1 + 2n), where n is an integer. 

Our main objective in what follows of this section is to analyse the cosmological 
evolution of the brane and, in particular, the singularities that may arise in a brane 
filled with DPGCG in the RSI model! 



3.1. DPGCG I 

We consider a brane filled with DPGCG I. In this model, the minimum energy 
density must be smaller than — 2A (see Fig. i.e. 2A < |^4| 1+2 ™, in order for the 
brane to have a Lorentzian evolution. This feature implies that the expansion of the 
brane starts at a non-vanishing scale factor defined^ as 



a\ 1 



1-(2A)TW|A|" 1 , (21) 



where the Hubble rate vanishes and p = — 2A. In the neighbourhood of oa i: the 
scale factor growth can be approximated by 

a * a Xl jl + 1(2X)t^\A\ [l - (2A)^ \A\^] (t - i Al ) 2 j , (22) 

which implies that the brane starts evolving smoothly, without any big bang sin- 
gularity. Moreover, the brane is born at a finite past cosmic time and a given 
radius a\ 1 , which is set up by the brane tension for a given DPGCG I, that should 
not be confused with a m i n given by Eq. ([3]), due to the fact that the brane tension 
acts as an early time energy cut off on the DPGCG I energy density. In other words, 
the brane tension excludes those radii of the brane such that a m j n < a < a\ t , where 
the energy density vanishes at a min . This cut off effect implies that the brane starts 
evolving at a finite past cosmic time as measured by an observer confined on the 
brane. 

Later on, the brane keeps expanding super- inflationarily (cf. Eq. (|17[) ) until the 
Hubble parameter approaches its maximum value Hi (or p acquires its minimum 
value; notice that 2A < \p\ < |^4| 1+2 ™), where 



f An anisotropic brane filled with a Chaplygin gas has been analysed in Ref|2Sl. 
g The analysis of this section is restricted to the expanding phase of the brane. 
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at large scale factor. It turns out that the brane is asymptotically de Sitter in 
the future despite that the DPGCG I energy density is negative and mimics a 
negative cosmological constant at large scale factor. At this respect, we recall that 
the modified Friedmann equation on the brane is quadratic on the total energy 
density of the brane (cf. Eq. (fT6|)) and therefore the effective energy density on the 
brane is positive . 




-2A 

1 

\A\T+ 

p 




Fig. 2. This plot shows the energy density p of DPGCG I versus the scale factor of the brane. 
The brane tension A acts as a cut off on the energy density at early time, therefore allowing only 
values of the scale factor larger than a\ ± , where p equates — 2A. At large values of the scale factor 
the DPGCG I energy density approaches a constant negative energy density for which DPGCG I 
mimics a negative cosmological constant, while the brane is asymptotically de Sitter. 



3.2. DPGCG II 

We move on to analyse what happens to a brane filled with a DPGCG II. First of 
all, let us point out that the main difference respect to the previous case is that 
the brane faces a big freeze singularity at a finite future cosmic time. Therefore, the 
fulfilment of the null energy condition 0^1 does not avoid the happening of such a 
singularity. 

In this model, the scale factor is bounded a\ 2 < a < a max , where a max is defined 
in Eq. (|T2")) and a\ 2 reads 



a\ 2 



l-(2A)Tw|^r 1 . (24) 



The brane starts its evolution at a non- vanishing value of the scale factor a\ 2 because 
the initial energy density of DPGCG II is larger than the energy scale fixed by the 
brane tension; i.e., | J 4| 1 /^ 1+a ' < 2A (see Fig. [3j). Therefore, A acts as an early time 
energy cut off on p and excludes the radii < a < a\ 2 where — 2A < p. Near a\ 2 , 
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the dependence of the scale factor on the brane proper time reads 

o(t) ~ a X2 jl + |(2A)t*^|^| [l - |^|- 1 (2A)^] (t - t A2 ) 2 | , (25) 

t\ 2 being the cosmic time at which the brane has its minimum radius. Unlike a brane 
filled with DPGCG III, where 2A < |^4| 1 ^ 1+Q \ here the brane tension prevents the 
existence of an infinite past for the brane. On the other hand, the Hubble rate grows 
so fast that the brane faces a big freeze singularity at a max . In fact, close to this 
doomsday, the cosmological expansion of the brane can be expressed as 

1- - \A\-—\n\-^ ' 1 (t max ~t)—i 



o(t) 



(26) 



where the big freeze takes place at i max and therefore both the energy density and 
pressure of DPGCG III blow up as well as the Hubble rate and its cosmic derivative 
do. The happening of this singularity can be understood as a consequence of the 
super-inflationary expansion of the brane (cf. Eq. (|17p ) and the unboundness of \p\ 
(2A < \p\). Moreover, the Hubble parameter grows at high energy; i.e. close to a max , 
and for a given constant a is enhanced by the brane effect respect to the standard 
4D situation represented by PGCG II and the PGCG IV fluids. Therefore, the big 
freeze singularity would take place in a more virulent way on the brane as H would 
diverges at a greater rate. 

3.3. DPGCG III 

The late time cosmology of a brane filled with DPGCG II or DPGCG III is the same 
and hence we refer the reader to see the previous subsection for a more detailed 
analysis. However, the early time evolution of the brane is quite different (see Fig. [3]). 
In fact, the brane tension does not exclude those radii banned by DPGCG II and the 
scale factor can take any value on < a < a max . The brane is then allowed to have 
very small radii because 2A < |^4| 1 /( 1 +") and therefore p < — 2A, a fact implying 
a well defined Hubble parameter even at very small scale factors (cf. Eq. (|16pV In 
other words, the energy scale set up by the brane tension is so small that it does 
not imply any cut off on p. It turns out that the brane evolution starts at an infinite 
time in the past. In fact, at small scale factor, the Friedmann equation implies 

a ~ a Q exp (Hi(t - t j) , (27) 

if 2A < |^4| 1 /( 1+Q ). In the previous expression Hi corresponds to the primordial 
Hubble parameter and is defined by Eq. (|23p and clq stands up for a small scale 
factor. Therefore, the brane is asymptotically de Sitter in the past. On the other 
hand, for the limiting case 2A = | J 4| 1 /( 1+Q \ the early time evolution is different (a 
sort of power law inflation) 



«o 




to) 



2(1+2") 
3 



(28) 
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a \-2 a max CL 




Fig. 3. This plot shows the energy density p of DPGCG II and DPGCG III versus the scale factor 

of the brane. The late time behaviour of a brane filled with any of these fluids is the same and 

shows a big freeze singularity. However, the early time evolution is different for the two fluids. For 

DPGCG III the brane tension A does not act as a cut off on the energy density (at early time) as 
i 

2A < | A| 1 + a . The reason is that the brane tension is too small. As the brane tension is increased 
(represented schematically by an arrow on the plot), the system behaves as a DPGCG II, where 
now A acts as a cut off on the energy density. In this case, the brane is allowed to have only values 
for the scale factor on the range < a < a. 



however, the brane has still an infinite past. We remind at this respect that 1 + In 
is negative. A summary of the results of the present section is given in table II. 



Table 2. This table summarises the asymptotic behaviour of a universe filled with each of the different types of DPGCG which were analysed in 
the text. 



A,B 


1 + a 


A 


a 


P 


Past 


Future 


A < 
B > 


(1 + 2n) _1 > 


2A < 


a \i < a < oo 


2A < \p\ < \A\W +a > 


finite past 


asymptotically de Sitter 


(l + 2rt)~ i < 


2A < \A\ 1 ^ 1+a > 


< a < a max 


\A\ l ^ x+a > < \p\ < oo 


asymptotically de Sitter 


big freeze singularity 


2A= \A\ L ^ 1+a > 


< a < a max 


|4|VU+«J < \ p \ < oo 


power law/infinite past 


big freeze singularity 


l^ji/U+aJ < 2A 


a^ 2 < a < a max 


2A < \p\ < oo 


finite past 


big freeze singularity 



4. The plain GCG with singularities 

One could think that the singularities appearing in section O in particular the big 
freeze one, could be a consequence from the use of a phantom fluid, since such a 
stuff violates the null, strong and weak energy conditional ^21 w e next show that 
even for GCG fulfilling such energy conditions a big freeze, a sudden or other type 



h A similar statement would apply to a dual PGCG as it does not satisfy the weak energy condition. 
Notice that a dual PGCG amazingly fulfils the strong energy condition. 
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of singularities may appear. This can be the case for the following set of GCG 
parameters: 

(i) A < 0, B > and 1 + a > 0. 

(ii) A < 0, B > and 1 + a < 0. 

(iii) A > 0, B < 0, 1 + a > and (1 + a) -1 = 2n, with n some positive integer 
number. 

(iv) A > 0, B < 0, 1 + a < and (1 + a) -1 = 2n, with n some negative integer 
number. 

On the other hand, we will also show that the dominant energy condition i s 
however not fulfilled at the big freeze singularity. We present a summary of the 
results found on this section in table III. 



4.1. GCG I 



We start analysing GCG I which is characterised by A < 0, B > and 1 + a > 
in Eqs. |T]) and @. This fluid behaves like a "dust" for small scale factor [J i.e. 
p/ p — > 0. Moreover, the universe starts evolving with an infinite Hubble parameter 
and an infinite deceleration parameter (a is very large and negative). Afterwards, the 
universe starts expanding, although never accelerating, until it reaches its maximum 
size given by Eq. (fT2|) . Now, how long does it take for such a universe to reach its 
maximum size? It can be easily checked that close to a max the scale factor growth 
can be approximated by 



1 - 



1 - 



«o 



l + 2a 
2(1 + a) 



3|A|(l + a-) 



for a ^ — 1/2 and by 



1 



1 



exp 



:H(*-*o) 



(29) 



(30) 



if a = —1/2. In the previous equations ao = a (to) is a scale factor close to a max . 

Consequently, if —1 < a < —1/2, the universe would reach its maximum size 
in an infinite cosmic time. Then the universe does not hit a future singularity. In 
fact, the scalar curvature is well defined and vanishes when the universe reaches its 
maximum size where the energy density and the pressure of GCG both vanish. 

However, for —1/2 < a the universe would reach its maximum size in a finite 
cosmic time in the future 



tr 



to 



2(1 + a) 
1 + 2a 



[3|A|(l + a)]- 



1 



no 



1 + 2q 
2(1 + = ) 



(31) 



1 Although, p/p — ► the pressure can diverge at small scale factors if —1 < a < 0. On the other 
hand, if a is positive, the pressure vanishes when the scale factor is small. 
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where it will face a curvature singularity if a is positive. Indeed, the scalar curvature 
blows up (see Eq. (E)) as well as it does the pressure of GCG. Therefore, there is 
a sudden singularity HI, a big brake singularity or a Type II singularity in the 
notation of Ref.^3 at t = £ max , with the peculiarity that the energy density vanishes 
at the singularity rather than approaching a positive value. If —1/2 < a < 0, the 
situation looks a bit more complicated. In this case, although the GCG energy 
density and the pressure are well defined at i max , and consequently also H and H 
do, it turns out that higher derivatives of the Hubble parameter can diverge when 
the universe reaches its maximum size. The n-th derivative of the Hubble parameter 
respect to the cosmic time has the structure given in Eq. ([9]). Therefore, following 
a reasoning similar to the one used in subsection 12.11 it can be straightforwardly 
shown that there always are higher derivatives of H that diverges at t = t max , unless 
a can be written as a = 1/ (2p) — 1/2, where p is a positive integer. This corresponds 
to a Type IV singularity in the notation of Ref. (please, see footnote |b]). In the 
particular situation where a — l/(2p) — 1/2, being p a positive integer, all the 
derivatives of H are well defined at t max and consequently there is no singularity 
at o max - 

Before concluding this subsection, we would like to point out that the universe 
never accelerates in this case. In fact, w = pj p is always positive. Therefore, the null, 
strong and weak energy conditions are always fulfilled. Nevertheless, the dominant 
energy condition is violated for scale factor larger than adomi where 



2 3 d+°)a n 



(32) 



consequently, the dominant energy condition is violated at the sudden singularity 
(see Fig. 0}. 





Fig. 4. Plot for the pressure/energy density ratio; w, for GCG I and GCG III, as a function of 
the scale factor. The dominant energy condition is not fulfilled for a scale factor larger than a;j om i . 
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Fig. 5. A schematic plot of a past big freeze singularity.This singularity takes places at a finite 
past cosmic time where not only the energy density but also the pressure blows up. Therefore, 
both the Hubble rate and its cosmic time diverge as well at the big freeze event. 



4.2. GCG II 

We briefly analyse what happens to a homogeneous and isotropic universe filled 
with a GCG II; i.e. with A < 0, B > and 1 + a < 0. It turns out that the 
scale factor is bounded from below, a m i n < a, where the minimum scale factor is 
given by Eq. ([3]). The GCG energy density as well as its pressure both diverge at 
amin, leading to a past curvature singularity (cf. Eq. ((8])). It corresponds to a big 
freeze singularity or a Type III singularity in the notation of Ref. (see Fig. [SJ. 
This singularity takes places at a finite past (see Eqs. (j4|) and (J6j) which also apply 
to in this situation). After the big freeze singularity, the energy density p keeps 
decreasing until it vanishes at very large value of the scale factor where the GCG 
has a dust-like behaviour. 

We would like to point out that, similar to the previous case, a universe filled with 
this fluid would never accelerate. In fact, w = p/p is always positive. Therefore, the 
null, strong and weak energy conditions are always fulfilled. However, the dominant 
energy condition is violated for scale factor smaller than <z<iom2, defined as 

«dom2 = 2" 3 < 1 +°)a min , (33) 

which implies that the dominant energy condition is violated at the big freeze sin- 
gularity (see Fig. [S]). 



4.3. GCG III and IV 

Finally, we analyse a universe filled with a GCG such that A > 0, B < and 
1 + a = l/(2n) where n is an integer. 

Let us consider that the fluid is characterised by a positive n which we will refer 
to as GCG III. Then, the scale factor is bounded from above by a max introduced in 
Eq. (fT2|) . At small scale factors the fluid behaves like "dust"; i.e. w ~ although 
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W 



1 








a, 



'mm 



a dom2 



a 



Fig. 6. Plot of the ratio of the pressure and the energy density; i.e. w, for GCG II and GCG IV, 
as a function of the scale factor. The dominant energy condition is not fulfilled for a scale factor 
smaller than ad om 2- 



p diverges, and therefore a big bang singularity shows up. Consequently, a FLRW 
universe filled with such a fluid would start at a curvature singularity. At a = a max , 
the energy density and the pressure of GCG III vanish and w acquires very large 
values, although the scalar curvature is well defined (see Eq. ([8])). Using Eqs. (|29|) - 
(|30| which hold, it can be proven that it will take an infinite cosmic time for a 
universe filled with this fluid to reach its maximum size. As it can be realised, 
the advantage of having a quantised a is to remove the possible different future 
singularities which may be induced by GCG I. After all, this advantage is again 
not that surprising as having n > implies —1 < a < —1/2, precisely a situation 
without future singularities for GCG I. Despite the absence of a future singularity at 
Omax and the fulfilment of the null, strong and weak energy conditions, the dominant 
energy condition is violated for scale factors larger than adomi (see Eq. (|32|) and 



For completeness, we briefly analyse what happens if n < a case which we 
name GCG IV. The scale factor is then bounded from below; i.e. a m i n < a, where 
flmin was defined by Eq. ([3]). At the minimum scale factor, p and p (as well as w) 
diverge inducing consequently an initial big freeze singularity (see Fig. [5|) because 
the singularity happens at a finite past cosmic time (Eqs. (j4]) and ((6]) both apply 
after replacing A{\ + a) by \A{1 + a)\). Afterwards, the universe starts expanding 
although never accelerating, and therefore whereas the null, strong and weak energy 
conditions are satisfied, the dominant energy condition is violated at radii smaller 
than adom2 defined by Eq. (|33| (see Fig. [6]) . On the other hand, at large scale factor 
GCG IV mimics a dust fluid. 



Fig.HD. 
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5. Discussion and conclusions 

In this work we have analysed in some detail a GCG-filled universe, pointing out the 
different kinds of singularities that may face in its future or past (see the tables for a 
summary of our results). In particular, we have shown that, although there are no big 
rip singularities in GCG models, other kinds of singularities can take place for other 
more complicated models. The emergence of some of this singularities, in particular 
the big freeze singularity, could be thought to be caused by the non-fulfilment of the 
null, strong and weak energy conditions by a phantom GCG0. For this reason, 
we have studied models of GCG satisfying the above mentioned energy conditions. 
We have then shown that even in this c ase a sudden ^\ a past big freeze and a type 
IV singularity, in the notation of Ref. could surface. Therefore, the fulfilment 
of the null, strong and weak energy conditions does not avoid the appearance of 
such singularities, particularly the big freeze one. We have also shown that in this 
case the dominant energy condition is not fulfilled at these sort of singularities (see 
Figs. Hand©. 

In section [3] we have shown how a brane can introduce a cut off for the energy 
density, for which case it modifies the model in such a way that the universe be- 
come always finite in the past (see Table II). As it is expected, brane effects can 
generally be relevant at sufficiently early times and, in fact, PGCG and DPGCG 
have late times behaviour which are rather indistinguishable, so that the universe 
essentially will become either asymptotically de Sitter or it will show a future big 
freeze singularity (see Table I and II). 

Before closing up, we should remind that close to a big rip singularity or a big 
brake even t, th e classical space-time breaks down and therefore a quantum analysis 
is required 1221. We similarly expect that quantum effects became also important as 
one approaches a big freeze singularity. In particular, the kind of past or future big 
freeze singularities we have point ed o ut in this paper need to be analysed from a 
quantum-mechanical point of view^S As a starting point, it would be interesting to 
perform such an analysis in the cases of the GCG II or GCG IV models, introduced 
in section|4j where at least we know that the null, strong and weak energy conditions 
are satisfied. We hope to report on this interesting issue in the near future. 
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Tabic 3. This tabic summarises the asymptotic behaviour of a universe filled with the four types of GCG analysed 
in the text. The future asymptotic behaviours labelled by (1) and (3) corresponds to —1 < a < —1/2 and < a, 
respectively. The future asymptotic behaviour labelled by (2) corresponds to —1/2 < a < 0, where a cannot be 
expressed as a = l/(2p) — 1/2, with p a positive integer. If —1/2 < a < and a can be expressed as a = l/(2p) — 1/2, 
with p a positive integer, there is no future singularity and the universe reaches its maximum radius at a finite future 
cosmic time. 



A,B 


1 + a 


a 


p 


Past 


Future 


A < 
B > 


positive 


< a < a max 


< p < oo 


dust-like 


(1) no singularity /infinite future 

(2) type IV singularity 

(3) sudden singularity 


negative 


(Imm < a < oo 


< p < oo 


big freeze singularity 


dust-like 


A > 
B < 


(2n)- L > 


< a < a max 


< p < oo 


dust-like 


no singularity/infinite future 


(2n)- L < 


Omin < a < OO 


< p < oo 


big freeze singularity 


dust-like 



